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Abstract. In this paper, we present a general derivation of a modified fluctuation- 
dissipation theorem (MFDT) valid near an arbitrary non-stationary state for a system 
obeying markovian dynamics. We show that the method to derive modified fluctuation- 
dissipation theorems near non-equilibrium stationary states used by J. Prost et al., 
PRL 103, 090601 (2009), is generalizablc to non-stationary states. This result follows 
from both standard linear response theory and from a transient fluctuation theorem, 
analogous to the Hatano-Sasa relation. We show that this modified fluctuation- 
dissipation theorem can be interpreted at the trajectory level using the notion of 
stochastic trajectory entropy, in a way which is similar to what has been done recently 
in the case of MFDT near non-equilibrium steady states (NESS). We illustrate this 
framework with two solvable examples: the first example corresponds to a brownian 
particle in an harmonic trap submitted to a quench of temperature and to a time- 
dependent stiffness. The second example is a classic model of coarsening systems, 
namely the ID Ising model with Glauber dynamics. 
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1. Introduction 

It is a general rule that as a system gets smaller its fluctuations increase. As a 
consequence, in small systems (like a colloidal particle or a biomolecule), thermodynamic 
quantities like work [TJ |2] or heat are only defined in a statistical sense. Exact relations 
between the statistical distributions of these thermodynamic quantities, known as 
fluctuations relations, have been obtained about a decade ago. Such ideas have lead 
to the emergence of a new field, concerned by the specificity of thermodynamics for 
small systems and which has been called stochastic thermodynamics. 

Fluctuation relations hold very generally for a large class of systems and arbitrarily 
far from equilibrium [U [3j SI El EJ [T] . They provide fresh ideas for revisiting an old but 
central issue of statistical physics, namely the origin of irreversibility of macroscopic 
systems. Furthermore, within the linear regime, these fluctuations relations lead to 
various new modified fluctuation-dissipation theorems (MFDT), which are interesting 
and valuable extensions of the classical fluctuation-dissipation theorem [HI EJ QUI E] . 

Recently, three main routes have emerged to construct such generalizations of the 
fluctuation-dissipation theorem: 

• In the first route opened by L.F. Cugliandolo et al. [12] and continued by E. 
Lippiello et al. [13] and G. Diezemann [H], the response function is written as a sum 
of a time derivative of the correlation function (similar to the equilibrium FDT), 
plus an additive function, called the asymmetry, which vanishes under equilibrium 
conditions. A physical interpretation for this asymmetry has been missing for many 
years, until recently, M. Baiesi et al. [15] propose to interpret it using a new concept 
called frenesy. This frenesy contains the time symmetric part of the non-equilibrium 
fluctuations. 

• In the second route opened by T. Speck and U. Seifert [16], the modifications of the 
fluctuation-dissipation theorem can be related to the so-called local velocity, which 
originates in the local currents present in the non-equilibrium situation. This route 
was further extended and generalized by R. Chetrite et al. p2] who also provided 
the MFDT with a Lagrangian frame interpretation [T8]. These ideas have then been 
confirmed experimentally using colloidal particles confined to circular trajectories 
[T9| [20] . In the end, it appears that the first and second routes are closely related 
and can be unified through the introduction of stochastic derivatives [21] . 

• In the third approach developed by J. Prost et al. [22] (see also R. Chetrite 
et al. p2]), the modified fluctuation theorem valid near non-equilibrium steady 
states (NESS) takes the standard equilibrium form except that it involves a new 
observable, function of the non-equilibrium steady state. This new observable is 
the system stochastic entropy, which here must be evaluated in the NESS [231 121] • 

Surprisingly, while the first and second approach provide ways to construct 
generalizations of the fluctuation-dissipation for cases where the system is initially in an 
arbitrary non-equilibrium state, the third approach seems limited to systems close to a 
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non-equilibrium stationary state. One objective of the present paper is to show that the 
third approach can too be extended to systems close to a general non-stationary state, 
in a way which is closely related to the first approach. In view of this, it appears that 
the three approaches provide closely related formulations to generalize the fluctuation- 
dissipation theorem to non-equilibrium situations. 

Such generalizations could potentially lead to a broad range of applications to non- 
stationary or driven systems, such as glasses, spin glasses, coarsening systems, granular 
media and dense colloidal systems, for which violations of the fluctuation-dissipation 
theorem have been an active topic for many years [121 E3 HH [25]. Another more recent 
but promising field of applications of these ideas concerns biological systems [261 EU EE] , 
where the application of MFDT could possibly lead to new methods to probe these 
complex media. 

In the next section, we show how to derive a MFDT classically, first using linear 
response theory, and then using more recent methods based on fluctuation relations. 
We then discuss an interpretation of the MFDT within stochastic thermodynamics, 
according to which the MFDT can be formulated in terms of a particular form 
of stochastic entropy. In the third section, we illustrate our framework with two 
pedagogical examples. In the first example, a Brownian particle placed in an harmonic 
potential is submitted to a quench of temperature and to a time- dependent stiffness, 
and in the second exemple, a ID Ising chain obeying Glauber dynamics is submitted to 
a quench of temperature and then probed with a magnetic field. 

2. Modified Fluctuation-dissipation theorem (MFDT) for general 
non-stationary states 

2.1. Stochastic modelling and definitions 

In the following, we derive a modified fluctuation-dissipation theorem (MFDT) for a 
system which is initially (at time 0) in a general non-stationary state. The evolution 
of the system at all times is assumed to follow a continuous-time Markovian dynamics 
of a pure jump type [29]. The transition rate to jump from a configuration c to a 
configuration d is denoted w t (c, c'), where the subscript t indicates that we allow for 
time-dependent rates. We denote p t (c) the probability to be in state c at time t. This 
quantity obeys the unperturbed master equation, 




(1) 



c' 



which can be written equivalently 




(2) 



c' 



in terms of the time- dependent markovian generator L t (c',c) defined by 




(3) 



C' 
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At time t — 0, an arbitrary but given time-dependent perturbation h t is applied to the 
system, and we denote by Pt(c, [h t ]) the probability to observe the system in the state 
c at a time t in the presence of this perturbation. The notation [h t ] emphasizes that 
the dependence is functional with respect to the perturbation. The evolution of the 
system for t > is controlled by the generator , which is defined similarly as in Eq. [3] 
provided the rates w t (c,c') are replaced by perturbed rates u^'(c, c'). This generator 
can be expanded to first order in [h t ] 

= L t + h t N t . (4) 

In the following, it does matter whether the unperturbed dynamics is autonomous or not. 
If this dynamics is non-autonomous, i.e. proceeds from the application of a protocol, we 
still denote the generator by L t without specifying this protocol explicitly. In this case, 
the application of the perturbation h t for t > can be seen as an additional protocol. 

Let us introduce 7r t (c, h) as the probability to observe the system in the state c at 
a time t > in the presence of a constant (time independent) perturbation h, which 
obeys 

^Vcfc)=5>(^/o^v,c). (5) 



A key object for the following discussion is ir t (c,h t ), which is constructed from 
7Tt(c, h) by replacing the time independent constant h by the value of the perturbation 
at time t, namely ht = h(t). In the particular case where the perturbed dynamics with 
a constant h is time independent (i.e. L\ = L h ), the subscript t in n t (c,h) may be 
dropped, and n(c, h t ) becomes the "accompanying" distribution introduced in Ref. [30] . 
We emphasize that 7r t (c, h t ) depends only on the perturbation at time t unlike Pt(c, [ht]) 
which depends functionally on the protocol history of the perturbation. The dynamics 
of 7Tt(c, h t ) is given by 

j t (n(c,h t )) = J2Mc\h t )L^(c\c) + h t ^^ (6) 

d 

= J>*(c'A) [L t ht (c',c) -5(c',c)htd ht ip t (c,h t )\ , (7) 

c' 

where in the second equation one has introduced ipt{c, h t ) = — \mr t (c, h t ) and an initial 
condition po( c ) — -Po( c ?0) = 7r (c, 0) has been assumed. 

The response function associated with a dynamic observable A t (c t , h t ) reads, for a 
perturbation applied at an earlier time t' > 0: 

where (..)[h t i represents an average with respect to the perturbed dynamics, and S/Shf 
is our notation for functional derivatives. We derive below a general formulation of a 
modified dissipation theorem for non- stationary states, which reads for t > t' > 

R(t,t') = --^(d h Mct',h)\ h ^ Mct,ht)). (9) 
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This relation qualifies for a modified fluctuation-dissipation because the response 
function is now expressed in terms of a correlation function of observables with respect to 
the unperturbed dynamics, denoted (..) [TO]. Note that when A t (c t ,h t ) = dh t ipt(c t , h t ), 
the modified fluctuation-dissipation theorem takes a more symmetric form, derived in 
Ref . [22] , under some specific conditions: these conditions are that the initial state must 
be stationary and that the dynamics followed by the system at constant perturbation 
h must be time independent (Lj 4 = L h ). Thus, as announced in the introduction, the 
present derivation extends the result developed in this reference to a more general initial 
condition and a more general dynamics. 



2.2. Derivation of MFD T from linear response theory 

This section shows how the modified fluctuation dissipation theorem of Eq. [9] follows 
from standard linear response theory. Starting from the master equation given above, 
one can generate a Dyson- type equation for the perturbed propagator P(ct\dt'), which 
is a fundamental result of linear response theory [9j [TOj [30] . This propagator P(ct\dt') 
represents the probability to find the system in the state c at time t given that it was in 
the state d at time t' according to the perturbed dynamics, while p(ct\dt') denotes the 
corresponding propagator for the unperturbed dynamics. When taken to first order in 
[h t ], the Dyson equation [31] for the propagator reads: 

P(ct\c 0) = p{ct\c 0) + f dt'h t '^2p(ct\dt')N t ,(d',d)p(d't'\c 0). (10) 

We then multiply this equation by an arbitrary observable, which we denote here A t (c) 
as a shorthand notation for A t (c,h t ). After integrating over the initial distribution 
Po( c o)> one obtains 

(A t (c t )) [ht] = (A t (ct)) + f dt'h t ,Y^A t {c)p{ct,dt')B t ,{d) (11) 

where we used the notation p(ct, dt') = p(ct\dt')p t >(d), and B is the operator such that 
B t ,(d) = Pt ,(dy 1 J2Nti(c",d) Pt ,(d / ). (12) 

c" 

The main point of introducing B is that it allows to write the response function in terms 
of a correlation function of two observables with respect to the unperturbed dynamics 
ED]: 

R(t,t') = (B t ,(c t ,)A t (c t )). (13) 

We can now use an expansion of the distribution ir t to first order in h, n t (c,h) = 
7T i°\ c ) + h-Kf~\c). It follows from the master equation that the zeroth order solution is 
Ttf\c) = pt{d), while the first order solution is 

= ^ 0) ^ , n(d,c)+4 1 \c')L t (d,c). (14) 
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Thus, the observable B defined above can be written 

B t (c) = ^(cy^Ntic',^^), (15) 

c' 



i -^n t 



^L t ( C ',cK (1) ( C ')). (16) 

After substituting this in Eq. [131 one obtains 



R(t,t') = J2Mc)p(ct\c't>) -J2Mc",c')4\c") 

c,c> \ c" 



(17) 



This form can be further transformed using the property that the unperturbed 
propagator p(ct\c't') satisfies the backward Kolmogorov equation [321 

d t ,p(ct\c"t') = -Y,LAc",c')p(ct\c't'), (18) 

c' 

so that in the end 

R ^ = [T,Mc)p(ct\c"t')4\c") ) , (19) 

\c,c" J 

which leads to Eq. [9] after using the relation n^\c) / nf\c) = <9hhi7r t /(c, h)\h->o- 

2.3. Derivation of MFDT from fluctuation relations 

For each path trajectory, we introduce the following functional 

y t = [ h T d h Mc T ,h T )dr. (20) 



Jo 

This functional y t has already appeared in Refs. [211 E3] but in a different form. In 
the appendix A, we explain the connections between the different formulations. The 
advantage of writing y t in the form of Eq. [201 besides its simplicity, is that the similarity 
with the functionals introduced by Jarzynski [1] and Hatano-Sasa [7] is then very 
apparent. 

In the same spirit as in the seminal works of Jarzynski and Hatano-Sasa, we 
consider below averages over trajectories with a weight y t . To perform such averages, 
we introduce the joint probability to be in the configuration c at time t with a value y 
for the quantity y t , Pj(c, y), which is defined by 

P t (c,y) = {s(c-c t )6(y-y t )) [fH] . (21) 

This quantity obeys the following master equation 

c' 

which can be solved through Laplace transform. Denoting Pt(c, 7) = f dyP t (c i y)e~ iy , 
we obtain 

— ol — = 2^ P *( C ' 7 ) L * ( c ' c ) ~ h ^ — ^ Ptfai). (23) 
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Thus, the equation satisfied by Pt(c, 7 = 1) is identical with the equation Eq. (j7]) satisfied 
by 7i t (c,h t ). Furthermore, since Po(c,y) = P (c)5(y), the two functions have identical 
initial conditions Po(c, 1) = f dy Po(c, y) exp(— y) = Po(c) = vr (c, 0). Therefore, these 
two functions must be identical, in other words: Pt(c, 1) = ir t (c, ht). Using the definition 
of the Laplace transform, it follows from this equality that 

ir t (c,h t ) = (6(c-c t )e- yt )) [ht] , (24) 

an equation which can be called a Feynman-Kac formula [2U [35] . 

By multiplying Eq. [24] by an arbitrary observable A t (c, h t ) and integrating over c, 
one obtains the following generalization of the Hatano-Sasa relation 

(A t (c t , h t )e- yt ) [ht] = j dc 7r t (c, h t )A t (c, h t ) = (A t (c t , h t )) n , (25) 

where in the last equality (..) 7ri denotes the average with respect to ir t (c,h t ). In the 
particular case that the initial condition is stationary and provided that the perturbed 
dynamics at constant perturbation h is time independent, the Hatano-Sasa relation [7] 
is recovered from Eq. [25] in the particular case of a constant observable A t = 1. 

We now consider a small variation with respect to the perturbation h t . From the 
definition of y t in Eq. [20], it follows that this quantity is small, at least of order one in 
ht, and thus e~ yt ~ 1 — yt- Therefore, we have 

(A t (h t )) nt ~ (A t (c t ,ht))[h t ] - ( / h T d hT ip T (c T ,h T )A t (c t ,h t )dT) , (26) 

WO / [ht] 

~ {At{ct,h t ))[ ht ]- drh T (d h ip T (c T ,h)\ h ^ A t (c t ,h t )) , (27) 
^0 

where in the last equation, we have approximated the derivative with respect to h T by 
a derivative with respect to h, an approximation which is correct to first order with 
respect to the perturbation, and at the same order in perturbation we have replaced 
the perturbed average by an unperturbed one. Taking into account that the functional 
derivative of the l.h.s. of Eq. [27] with respect to h t > vanishes for t' < t, and rewriting 
the second term of the r.h.s using an integration by parts, we obtain Eq. [9] 

2-4- Interpretation of MFD T in terms of trajectory entropy 

In section |2~T| we have introduced a key quantity namely ipt{c, h) = — In 717 (c, h). When 
properly evaluated on a specific trajectory [c t ,h t }, the function ip t (c t ,h t ) gives access 
to the functional 3^ defined in Eq. [2Q1 and when evaluated at small constant h on a 
trajectory [c t ], it allows to calculate the response in the MFDT according to Eq. [9] 
Clearly, this quantity must be closely related to the stochastic entropy introduced in 
Ref. [36]. Indeed, the stochastic system entropy is defined as 

8t(c t ,[h t ]) = -hxPt{c t ,[h t ]), (28) 

and therefore depends functionally on the perturbation [ht]. In contrast to that, the 
system entropy which is needed here is evaluated using a constant perturbation [h] , 

s t (c t , [h]) = -\nP t (c t , [h]) = -hxTT t (ct,h) = Mct,h). (29) 
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We now focus on the trajectories taken by the system, which can be described by a 
set of discrete values C = {c , ci..cjv}, with the convention that the system is in state c 
at time and in state cat at time t. Furthermore, the transition from state Cj_i to state 
Cj occurs at the jumping times Tj. The stochastic system entropy can be decomposed 
as St(ct, [h t ]) = — s[(c f , [h t ]) + [h t ]), in terms of the reservoir entropy production 

[h t ]) (also called medium entropy in Ref. [36]) and the total entropy production 

4 0t (c t ,[h t }). 

The system entropy is a state function, which means that 

As t (c t , [ht]) = -\nP t (c u [h t ]) +mP (coA)- (30) 

In contrast to that, the reservoir entropy and the total entropy are not state functions, 
but are trajectory dependent quantities, which can be written 

AsXct, [ht]) =£ln^^, (31) 

j=l PTjiCjAhtUWr-iCjiCj-x) Jo 

where hj is the value of the control parameter at the jump time Tj. Between the jumps, 
is a constant function of the time while s t and s* ot are in general non-constant but 
continuous functions of the time. All these functions, St, s r t and sl ot are discontinuous 
at the jump times Tj. 

When adapted to the case of a constant perturbation [h] , the above decomposition 
of the system entropy leads to two terms in the MFDT. Starting from Eq. [9] together 
with Eq. |30l one obtains 

R(t,t') = (d h Astict>, [h])\ h ^ At(ct)) = R eg (t,t') - R neq (t,t'), (32) 
R eq (t,t') = A (d h As r t ,(ct>, [h])\ h ^A t {ct)) , (33) 

Rne q (t,t') = -^{d h A S y(ct,, [h])\ h ^ A t (ct)) . (34) 

This decomposition of the MFDT contains two terms, the first term R eq (t,t') which is 
analogous to the equilibrium FDT in a sense made more precise below, and the second 
term R neq (t,t') which represents an additive correction. Such a decomposition has been 
discussed by several authors following the original work of Ref. [IB]. Note that the 
interpretation of the MFDT or standard FDT in terms of trajectory entropies is more 
recent [151 EZ]; in the context of non-equilibrium stationary states this has been done 
in [23| [23] . The present decomposition is very similar to that case, but here the total 
entropy production replaces the so-called adiabatic entropy production [38] , because the 
non-adiabatic part is non-zero and contributes to the second term R neq (t, t'). 



where 



and 
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Reformulation of the MFDT using local currents Let us rewrite more explicitly the 
two terms above of the MFDT without the time derivatives present in Eqs. I32H341 but 
instead using local currents. For the first term, using the expression for the reservoir 
entropy of Eq. [3TJ we have 



rt' 

(As r t ,(c t ,,[h])A t (c t ))= / ds Ps(c)w s (c,c')ln 

^° c,c',c" 



v w h s {c,J) 



which implies that 



d 



w%(c', c) 

-jp {dh As[,(c t /, [h])\ h ^ A(c t )} = ^ft-fcKtc.c'Rlnl^ 

c,c',c" W t'V C i C ) 

= {jvMMct)) = Re q (t,t'), 

where jt' denotes the local current given by 



x P (c"t\c's)A t (c"), (35) 



Jt 



Pt'(c) , ,s „ , u>iHc, c') 



p(c"t\c , t')A t (c"), 
(36) 

(37) 



h->0 



A property of this local current is that its average represents a physical current: 



\h->0 



(38 



with J t i(c',c) = p t >{c')w t >{d ', c) — /O t /(c)wt'(c, c'), the unperturbed probability current 
between the states c and c'. 

The same strategy can be used to transform the second term in the MFDT. 

v K s (c,h)w h 8 {c,c>) 



'M?(c t ,, 



[h])A t {ct)) = / rfs V p s (c)u7,(c,c') In 

+ / dr^p r (c')(^ r )(c',/i)p( C ' / t|cV)A( C ' / ) 



p(c"t| C ' S )A(c") 
(39) 



which implies 
d 
dt 



d Irs a M, r, n\l w \\ V^/ Pt'(c) / ,„ , 71*/ (c, k)wb (c, C+l) 

+ (^#)(c f /,/i)|^ A t (ct)) . 



Met)) 

(40) 



One can rewrite this in a more compact form in terms of another local current v t > such 
that 



d 



— (d h AafCc, [/i])| /w0 A(c t )> = M<*)^t(c*)> = R neq (t,t'), (41) 



where 



„ Pt'( c ) 



(42) 



Note that this equation together with Eq. [38l imply that both currents vt>(d) and jV(Q') 
have the same average. In the end, the MFDT takes the form : 

R(t,t') = (UtrM - M<*))M<h)) ■ ( 43 ) 
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We emphasize that the decomposition of Eqs. I3"7H4"21 is a general result, which does not 
rely on any assumption about the form of the rates. Note that the response function 
can thus be written in a way which does not contain a time derivative provided two 
local currents j t > and v t i are introduced. This comes at the price that there is not a 
unique decomposition of this type [24] . 

When additional assumptions are available about the transition rates, the two terms 
can be further transformed. One such assumption is a generalized detailed balance 
relation of the form : 

w f 'fc, c?) wAc,c') „ ... . . 

Wt(c',c) w t {c',c) 

where d t (c,c') represents the variation of a physical quantity conjugate to h t during the 
transition from state c to d at time t. In this case, the local current j t i can be simplified 

as : 

jAc') = ^ £ ^TJS^(c,c')d t (c,c'). (45) 

One particular usual choice of transition rates compatible with Eq. HH is 

w h / (c, c') = w t , (c, c>) exp Uhv \ , (46 ) 

where 0(c) represents a physical time independent observable, and d(c, c') = 0(c')—0(c) 
its variation between state c and d . See [2TJ for a general discussion on the various forms 
of perturbations in connection with the potential theory. 

The reservoir entropy introduced in Eq. ED only needs to be evaluated for a constant 
perturbation [h], therefore 

N h ( \ N 

d h Asl(c u [h])\ h ^ = T, d ^ W H( \ =/3 E (0(^-0(0^)), 

j=1 w n\ c h c 3-±) h^o j=i 

= p (0(c N ) - O(c )) = P (0(c t ) - O(c )) . (47) 
Substituting this in the first term of the MFDT, namely, Eq. [331 one obtains 

Req(t,t') = /3-^ (0(ct>)A(Ct)) , (48) 

which is a form similar to the equilibrium FDT. Indeed, in the case of a perturbation 
around an equilibrium state, the average (..) becomes an equilibrium average, and the 
equilibrium form of FDT is recovered. We explain in the appendix A how the second 
term R neq (t,t') can be transformed in the continuous space limit for the particular case 
of a nearest-neighbour random walk. 

In the end, the present framework provides a way to interpret the MFDT at the 
level of trajectory entropies in the case that a system is perturbed near an arbitrary 
non-equilibrium state. Thus, this framework generalizes the results obtained in Ref. [23] 
for non-equilibrium stationary states. 



11 



3. Illustrative examples 

While the derivations above concern discontinuous pure jump Markov processes, the 
results are more general and their proofs are transposable for continuous diffusion 
processes [TO] . For this reason, we provide in the following two illustrative examples 
of each kind: the first one corresponds to a continuous process of the Langevin type 
while the second one corresponds to a discrete jump process (Glauber dynamics). 



3.1. A particle obeying Langevin dynamics in an harmonic potential with time 
dependent stiffness and submitted to a quench of temperature 

We consider a particle in one dimension and in an harmonic potential obeying Langevin 
dynamics: 

7 7 ^ r J~' 

x t = — -x t + — + r] t with (r) t r) t >) = — -5(t-t') and (r) t ) = 0, (49) 
7 7 7 

where rjt is a Gaussian white noise, k t a time dependent spring constant, 7 a friction 

coefficient and T t the time dependent temperature of the bath, which starts from To 

at t — and ends at T t at t — tf. As a result of the non-stationary bath and of the 

time dependent spring constant, the system at time t > is not in equilibrium although 

it is assumed to be at equilibrium at t = 0. We denote by h t an additional external 

perturbing force. For this system, one can compute explicitly, provided that the spring 

constant is integrable on interval [0,t], the position at time t, 



x t = x e~ drK ^ + e ~ $ drKh / dr 



1 



Jl dT 'K<h 



(50) 



which is also a random Gaussian variable thanks to the linearity in 77 and because 
the probability distribution of x is the equilibrium Gaussian one at T , i.e po(x ) = 
exp(— koXft/ (2Tq)) / Z) with Z the partition function. From Eq. [501 we obtain 

dr — exp 



2 



a 



( x t)[h t ] = 
( x t)[h t ] ~ 
(xl) exp 



7 



[hi 



drkr/'j 



dr'k T i 1 7 



2T T 
dr — - exp 

7 



dukuj^ 



(51) 
(52) 
(53) 



The functional derivative of Eq.[5T]with respect to h t > gives directly the response function 



R{t,t') 



8{x t ) 



[hi 



5h+i 



exp 



drk T j^i 



(54) 



h^o 7 

As we show below, this result can also be recovered from the MFDT of Eq. [9j Since Xt 
is Gaussian variable, we deduce from this that the probability density function at time 
t is 



Pt(x, [ht]) 



1 



exp 



-1 

2oi 



x 







dr — exp 

7 



t 



(55) 
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We note that of does not depend on the perturbation [h] but fi t does. To obtain 
the probability density function n t (x, h), we just make the perturbation constant using 
ir t (x, h) = P t (x, [h]), so: 

ft / rt n n 2 " 



2of 



^l dTexp {-L dT%,/i )) 



(56) 



Note that 7r 4 (x, /it) is indeed distinct from Pt(x, [ht]) as emphasized from the beginning. 
Since we now have the key function ir t , we can compute the response function and the 
functional y t . 

To obtain the response, we calculate ipt{x, h) = — lmr t (x, h) and its derivative with 
respect to h: 

d h ip t (x, h) = — \{x - hl t )l t , (57) 
of 



where 



resulting in 



I t = - /Ve(-/>' fc ^), (5? 



7 Jo 



(d h ^t>{x t >)x t )\ h _, Q = -% (x t x t >)\ h ^ ■ (59) 



This last expression requires the two times correlation function for x that can be obtained 
from Eq. \5U\ assuming t' < t 

(x t >x t ) = HtHv + of, exp (- drk T /j\ . (60) 

Since that fi t = fit' = in the limit h — > 0, we obtain 

(d h ip t '(xt')xt}\ h ^ = -If/expi- / drk T /^f) , (61) 



-t' 



t> 



dTe( f ° dT ' k T'/"<)exp(- drk T /^j . (62) 



Uo 

Now, by taking the opposite of the time derivative with respect to t', we obtain the 
response function of Eq. H3J which provides a verification of Eq. |9] on this particular 
example. Note that the response function does not have time translational symmetry 
(it is not solely a function of the time difference t — t') in the general case that the spring 
constant is time- dependent. Furthermore, this response function has the property of 
being independent of the protocol of temperature variation. 
We obtain y t from Eq. [20] and Eq. E3 

t Kir 

y t = - / dr- L ^-{x T - h T I T ), (63) 

Jo <*i 

which is linear in x T and hence is also a Gaussian variable. Its mean value and variance 
of y t are respectively 

(y t ) = fdr^iKIr-^), (64) 

./n 0_ 
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Figure 1. Probability density functions Pt(y) at five different times t for a particle 
in an harmonic trap obeying Langevin dynamics. The symbols represent an estimate 
of Pt(y) based on 10 5 trajectories of total duration tf = 5.12; the solid line is the 
Gaussian probability density which has the mean and variance given respectively by 
Eq. [64] and 1651 The system is at t > in a non-equilibrium state due to an imposed 
time-dependent spring constant k t — 5 + 2.5sin(|j), and a time-dependent heat bath 
temperature, which is such that it is To = 5 at t = and T t = 1 for t > 0. The 
friction coefficient is 7 — 1. This system is further perturbed by a force, according to 
the protocol ht = 5sin(|^). 
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Figure 2. Numerical test of the detailed fluctuation relation of Eq. [66] satisfied by the 
probability distribution of y. Here t = tf = 5.12, and all the other parameters are the 
same as in Figure [U 

{y 2 t ) - (y t ) 2 = f dr f dr' ^p^ ((x T x T ,) - ^ r /v) • (65) 
Jo Jo a T a T ' 
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In the end, after transforming Eq. [65] using an integration by parts, we find that 
(y^) — (yt) 2 = 2 (3^t) • Since y t is Gaussian process, this relation implies that the 
probability density function Pt{y) satisfies the following detailed fluctuation relation 

p~(~y) =exp{y) ' (66) 

which we have also confirmed through a numerical determination of the distribution of 
y t as shown in figures [TH21 

Remarkably, this detailed fluctuation relation holds although no dual process has 
been invoked. Thus, this relation differs in an essential way from the generalized 
Crooks theorem given in [21] (see Eq. 133). Another difference between Eq. [66] and the 
generalized Crooks theorem of that reference is that in the generalized Crooks relation, 
the initial condition in the forward or backward processes are taken according to the 
distribution ttq(x, ho) and TTt(x, ht) respectively. In contrast to this, in Eq. [66] the initial 
condition in which the system is prepared results from the application of unperturbed 
dynamics. For this reason, one may say that Eq. [66] is closer to a relation of the Bochkov 
and Kuzovlev type rather than to a Crooks relation [39J. 

Finally, we would like to emphasize that Eq. [66] is a very general result for linear 
Langevin dynamics. We have checked in Appendix B that Eq. [66]can indeed be extended 
to a general multidimensional linear Langevin dynamics. 



3.2. The ID Ising model with Glauber dynamics 

3.2.1. Introduction We now move to a more complex system with many interacting 
degrees of freedom, which will allow for phase transitions and ordering phenomena 
absent from the previous example. The system is the Ising-Glauber chain in ID, which 
when submitted to a temperature quench, is a paradigm for coarsening dynamics [4"0] . 
For this system, explicit exact expressions of the correlation and response functions have 
been obtained; and the ratio between these two quantities admits a non-trivial limit, 
which is a universal quantity in the case of a quench to the critical temperature |41j . 
Multi point correlation functions have also been calculated analytically in order to test 
theoretical ideas about the dynamic heterogeneities of glasses [121 [33]. In more general 
spin systems, the correlations or response can not be obtained analytically, but the 
response function has been shown to be related to correlations characteristic of the non- 
perturbed system [131 CLU El US] , a conceptual progress but also a definite advantage for 
numerical simulations as compared to previous methods. In the following, we illustrate 
the framework of modified fluctuation-dissipation theorem presented in the previous 
sections for the Glauber-Ising chain submitted to a quench of temperature. Using 
analytical calculations, we first show that we can recover the known exact response 
function using this formalism. We then present some numerical simulations to confirm 
the theoretical expectations. 
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3.2.2. Definition of the rates This Ising-Glauber chain is made of L Ising spins cr^ = ±1 
with i = 1..L in one dimension, and is described by the following Hamiltonian 

L 

^(W) = - J 5^ ^^+1 ~~ H rnO- m , (67) 

1=1 

where J is the coupling constant and H rn a magnetic field which acts on the spin m. We 
assume periodic boundary conditions. The magnetic field H m will be the only control 
parameter. In principle, we could allow for many control parameters corresponding to 
magnetic fields present on any lattice site, but since we are mainly interested in the 
linear response regime, we restrict ourselves to the case where this magnetic field only 
acts on the spin m. The Ising chain is assumed to be initially in equilibrium at an 
infinite temperature. At t = 0, it is submitted to an instantaneous quench which brings 
the temperature to T, and the system in a non-equilibrium state. The system further 
evolves after the quench from the time to the time t' where the small magnetic field 
H m is turned on. Therefore we assume that H m {t) = H m 9(t — t'), with t > and 6 the 
Heaviside function. 

The probability to find the system in the state {a} = {o~i, <jl} at time t, Pt({a}), 
obeys the following master equation 

i i 

where w Hm ({o~}, {c} 4 ) is the rate to jump from the configuration {a} to the configuration 

WY = {°i> -,o-l}- 

Following Ref. [JT], we choose the rates in the presence or absence of a field to be 
given respectively by 

w Hm ({a}, {a} 1 ) = |(1 - a i tanh(/3J(a i _ 1 + a i+1 ) + (3H m 5 tm )), (69) 

MM, WY) = f (i - ^i-i + , (70) 

where a denotes the inverse characteristic time scale of the transitions (which we take 
below to be equal to 1), 7 = tanh(2/3J) and (3 = 1/T, the inverse of the temperature 
after the quench. Note that these rates depend on time (for t > 0) only via the control 
parameter H m . It is important to also realize that this form of the rates is just one of 
the possible choices compatible with the detailed balance condition, which imposes that 

"> fl (W,K}) _ e- W) m , 
w H ({a'},{a}) e-WW) ' 1 1 

Other forms are possible, for instance Eq. H6] corresponds to a different acceptable choice, 
which is the one made in Ref. [13] while Eqs. [oWTUl is the choice of Refs. [HI US] . 

3.2.3. Analytical verification of the MFDT Unlike in the previous example of a particle 
obeying Langevin dynamics, where the non-stationary probability distribution 7r t was 
analytically solvable even in the presence of the perturbation due to the assumption of 
an harmonic potential, in the present problem, the probability distribution 7r t ({cr}, H m ) 
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to find the system in an arbitrary configuration {a} at time t with a constant magnetic 
field H m applied from t = is difficult to obtain even in the absence of a magnetic field 
(H m = 0). As we show below, this distribution is not required to evaluate the response 
function for a general one spin dependent observable A(a n ), because in this particular 
case only the reduced one-spin distribution 7r t (<7 n , H m ) = J2{ a }=ia n ^(l "}' Hm) is needed 
and fortunately, this reduced distribution can be calculated analytically We define the 
response function as 

5(A(a n (t))) [Hr , 



Rn—m(ty t ) T 



(72) 

ffm^D 



6H m (t>) 

which contains an extra factor T with respect to the definition used in previous sections. 
The reason for this extra factor is purely a matter of convenience, but historically it 
was introduced in order to provide the response function with a well defined limit when 
T — > 0. According to the MFDT of Eq. [HI this response function is 



dP 



T^J2 A MpKt\{a'}t')d Hm MW},H m )\ 



0",i,{o"'} 

where in the last step, we used 7iv({er'},0) = pv ({&'})■ 

To progress, we need an explicit expression for the propagator p(a n t\{a'}t') , which 
is the probability to find the system with spin n in state a n at time t in the unperturbed 
dynamics given that the system was in the state {a'} at time t'. In fact, this propagator 
is directly related to the average magnetization at time t when the system starts in the 
state {a 1 } at time t', namely (cr n (t)){ cr }( t /) = { cr /}, 

p(a n t\{a'}t') = - ((T n ((Tn(t)){a}(t')={a>} + l) • (73) 

From Eq. [68l one can show that this average magnetization in the absence of an applied 
field satisfies jlQ] 

(v n (t)){a}(t')={*>} = ^2 G ^-k(t - t')a' k , (74) 

k 

where Gk{t) is the Green function of the problem. In the following, we consider the 
thermodynamic limit L — > oo, in which case, the Green function can be written as 
Gk{t) = e~*/fc(7t), in terms of Ik the modified Bessel function. It follows from the above 
two equations that 

p(a n t\{a'}t') = ^J2°'k G n-k(t-t')+ 1 -. (75) 
k 



After substituting this expression in Eq. [72j one obtains 

dt 



Rn- m {t,t') = T^- 5>K) d Hm 7r t/ ({a'},H m )\ Hrn ^ 



(76) 
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The fact that the observable A(a n ) only depends on a single spin a n leads to a 
simplification of the expression above since it is possible to sum over all spins in {a'} 
except for the kth spin of that set. The second term in the bracket of Eq. 1761 vanishes 
due to the normalization condition Ylfa'} ^'(W}, H m ) = 1; and we obtain 

R n _ m (t,t')=T^Y. E \v' k VnA{a n )G n - k {t-t') d Hm Tr t ,(a' k , H m )\ Hm ^ .(77) 



k <J\,<J n 



This equation shows that the response now only depends on the reduced one-spin 
distribution 7i t (a n ,H m ) = J2{a}=£a n 7l t({o'},H m ), which can be obtained analytically at 
first order in the applied field H m from the magnetization. Indeed, in the presence of a 
field, the magnetization at a time t > is 

M*)> W(o)=K} - Y, G»-fc(*K + flH m [ G n . m (t - t')K(t')dt\ (78) 

Jo 



k 

2 



7 

where K{t) = 1 — (1 + (a n+ x(t)(T n -i{t))), a correlation function which is known 

analytically. For instance, in the case of a T = quench where 7(t) = 1 for t > 0, we 
have K{t) = e- 2t (I (2t) + I 2 (2t) + 2/ 1 (2t))/2 06]. 

Since the chain is initially at infinite temperature, the average magnetization of 
spin k at time t — vanishes for all k. It follows that the first term of Eq. [751 vanishes 
when the average over the initial condition is performed. Using Eq. [75l we obtain the 
distribution 7i t (a n , H m ) at first order in the applied field 

7r t {a n , H m )~^ + ^(3H m a n J G n ^ m (t - t')K(t')dt'. (79) 

We can use this distribution to write the response function in Eq. [77J in a more explicit 
form. After summing over a' k , we have 

d 1 /"*' 
R n -m(t, t') = — ^ -a n A(a n )G n _ k (t - t') J G k ^ m {t' - u)K(u)du 



A(l)-A(-1) d 



2 dt 
A{l)-A{-\ 



r? 

- / V G n _ fc (t - H)G h ^ - u)K{u)du 
Jo , 



G n . m (t-t')K(t'), (80) 

which agrees indeed with the response function obtained in [JT] in the case A(a n ) = a n . 

Besides recovering the known response function of the Glauber-Ising chain, we can 
also investigate the separate contributions of the two local currents j t > and v t > introduced 
in the section 12.41 It is straightforward to show using the detailed balance condition 
of Eq. [71] that the first term in the MFDT, namely R eq (t,t'), can be written in a form 
similar to that of the equilibrium FDT (note that in the equation below, the index j 
represents the discrete times when jumps occur): 



d Hm Asl({a} t ,[h])\ Hm=0 = > d Hm ln 



<1W(i-i),W(i)) 



^ m ^< m (W(i) ; W(i-i)) 



(81) 

H m =0 
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/3^K(J)-<UJ-1)), (82) 



13 (a m (N) - a m (0)) = (3 {a m {t) - a m {0)) . (83) 



Substituting this in Eq. [33J one obtains the equilibrium part of the response for a general 
multi-spins observable A({a} t ), 



Alternately, one can also recover this result through an evaluation of the local current 
jt> using Eq. [371 

For the term associated with the local current v(t'), one can show either from the 
decomposition of the stochastic entropy into two terms or from Eqs. I4TH421 that the 
expected part of the MFDT is recovered, in other words that R neq (t,t') = R(t,t') — 
R e g(t, t'), where R(t,t') is the response given in Eq. [80] for the case of a one-spin 
observable A{a n ). In the literature [Ml [2], the term R neq (t,t') is called the asymmetry, 
its precise form depends on the specific form of the rates (unlike the first term which only 
depends on the ratio of forward to backward rates) and it vanishes under equilibrium 
conditions. The present derivation makes also clear that there is in principle a physical 
observable associated with this term, namely As t t ot ({a} t , [h]). 

3.3. Numerical verification 

As mentioned above, the distribution ir t ({a},H m ) does not seem to be accessible 
analytically. In order to test our framework, we have thus determined these distributions 
numerically from simulations for various values of H m and {a}. Then, we have calculated 
the response function R(t, t') via Eq. [9] and using trajectories which were simulated 
according to the dynamics in the absence of a magnetic field. In figures [3J1H the 
integrated response defined by 



is shown, where the symbols represent the response function obtained from such 
simulations at zero field via Eq. [9] and the solid line is the analytical expression obtained 
from Eq. [801 Since Eq. [80] takes a simple form in Laplace space [H], this solid line was 
obtained through a numerical inverse Laplace transform of that equation. 

In the simulations at zero field, we have used a small system size of L = 14. 
It is difficult to go to significantly larger sizes with the present algorithm, because the 
numerical determination of the distribution n t ({a}, H m ) becomes rapidly a difficult task 
in large systems given that the configuration space grows as 2^. 

4. Conclusion 



Re q (t,t') = P—(A({a} t )a m (t'))- 



(84) 




(85) 



In this paper, we have presented a framework which can be used to generalize 
the fluctuation-dissipation theorem to non-equilibrium systems obeying markovian 
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t 



Figure 3. Integrated response functions %(f, t') versus time t for a step protocol of 
magnetic field H m = 0.05 starting at various values of the waiting time t' after the 
initial quench at t = 0. The different values of t' arc t' = for circles, t' = 1.43 for 
squares, t' = 7.14 for triangles and t' — 14.3 for inverted triangles. The parameters 
are the following: J = 0.5, a = 1, T = 1 and L = 14. The response is calculated on 
the same spin where the magnetic field is applied, here n = m = 3. The averages have 
been done with 10 6 trajectories of 400 time steps of length dt = 0.07. The continuous 
lines stand for the integrated response obtained analytically, while the symbols have 
been obtained from the MFDT of Eq. [9] using the unperturbed dynamics. 



0.05 




Figure 4. Integrated response functions %(f, t') versus time t in the same conditions 
as in figure [3] except that now the response is evaluated on a different spin (n = 7) as 
compared to where the magnetic field is applied (n = 3) and the averages have been 
done with 10 7 trajectories. 

dynamics. We have shown, using alternatively the linear response theory or a first 
order development of a fluctuation relation, that the main result of Ref. [22] for systems 
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in a non-equilibrium stationary state is generalizable to systems which are near non- 
stationary states. This generalization is important because it restores a fluctuation 
dissipation theorem in the form of a unique correlation between physical observables 
and it opens many new possibilities to apply this framework to experiments. 

In fact, this framework is applicable to systems which are prepared in a non- 
equilibrium state, and which are then be further probed through the application of 
a time- dependent control parameter. This situation is typically the one encountered in 
studies of slow relaxing or aging systems, but it is also a frequent situation in biological 
systems. One outcome of our approach is that it is possible to replace this complicated 
problem by a somewhat simpler problem, namely the problem of determining the 
probability distribution n t (c,h) to find the system in a non-equilibrium state c but 
with a time-independent perturbation h. Having to consider only a time-independent 
perturbation to probe a non-equilibrium system, should be a definite advantage both 
from an experimental and theoretical point of view. 

Our study of a Brownian particle in an harmonic potential and submitted to a 
quench of temperature, raises the question of the validity of fluctuation relations for 
particles in contact with a nonequilibrium bath. This kind of studies may be important 
to understand the nonequilibrium fluctuations of a Brownian particle confined in a gel as 
in the experiment of Ref . jl?] . We hope that our work will trigger further studies on the 
applications of stochastic thermodynamics to the characterization of non-equilibrium 
systems, and in particular for non-equilibrium systems which result from contact with 
a nonequilibrium bath or from a coarsening process. 
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Appendix A. Link with previous formulations 

Equivalence between different expressions of the work like functional 

As mentioned in the main text, the functional [Vt defined in Eq. [20] has appeared before 
in Ref. [33]. We explain here how to make contact with the different notations. The 
functional denoted here 3^, corresponds to the one called Wt in Eq. 24 of this reference. 
In order to see this correspondence, one should choose the quantity defined as f t in 
Ref. |33j to be equal to n t (c, h t ). With this choice, 




dr I 7T r (c T , h T ) 1 (<9 r 7r r ) (c r , h T ) — (d T hi7r r ) (c T , h T 



d_ 




) - h T d hT lmi T (c T ,h T ) 



(A.2) 




) 
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where in the first line, the dagger stands for adjoint operator and the • indicates that 
the expression is to be understood as matrix product before evaluation on trajectory c T . 
This calculation shows that Wt indeed coincides with y t defined in Eq. 1201 

Equivalence between the discrete and continuous formulation of the MFDT 

In Refs. [171 [33], a modified fluctuation-dissipation has been derived for general 
continuous diffusion processes. Here we show that the present framework formulated 
for discrete jump processes leads to the same results, when the appropriate continuous 
limit of the master equation is taken [4"5] . 

We shall assume the same parametrization of the rates as that given in Eq. H6j The 
equilibrium contribution in the response, R eq (t,t'), can be written in the same way as 
in Ref. [33], therefore, we shall focus on the other non-equilibrium contribution namely 
R neq (t,t'). For simplicity, let us consider a nearest-neighbour random walk on a ID 
lattice, in which a random walker at position m can only jump to neighbouring sites 
mil. We denote the actual distance between all the sites by e. The master equation 
is 

^Pt(m) = J t (m+ l,m) - J t (m - l,m) (A.3) 

with J t (m,m + 1) = p t (m)w t (m,m + 1) — pt{m + l)w t (m + l,m). Using the Taylor 
expansion p t (m ± 1) = pt(m) ± d m pt(m), we can rewrite the discrete currents as 

Jt(m, m± 1) = Wt(m,m±l)p t (m) — iy t (m±l, m)(pt(m) ±d m pt(m)).(AA) 

To establish a link between these discrete currents and the current arising in the 
corresponding Fokker-Planck equation, we introduce the notations u t {m + 1) = 
w t (m,m + 1) — w t {jn + l,m) and 2D t {m + 1) = w t (m,m + 1) + w t (m + l,m), so 
that the discrete currents are 

J t (m,m + 1) = pt(jn)(u t (m) + d m u t (m)) - w t (m + l,m)d m p t {m) (A. 5) 
J t (m,m - 1) = - p t {m)u t {m) +w t (m- l,m)d m p t (m) (A. 6) 

In the continuous limit, e — ^ 

pt{m) ~ ep t (x), d m ~ ed x , u t (m) ~ u t (x)/e, D t (m) ~ D t (x)/e 2 , (A.7) 

and the discrete currents can be related to the current Jt(x) entering the Fokker-Planck 
equation: 

J t (m, m ± 1) ~ ±u t (x)p t (x) =F D t d x p t (x) = ±J t (x). (A. 8) 

Let us recall the expression given in Eq. [42] for the local current v t i(m\. 

J t ,(m,m+l)d h w$(m,m + l)\ h -> . . 

u t i{m) = A.9 

Pt'{m) 

J t ,(m,m- l)d h w$(m,m- l)\ h ^o 
Pt'{m) 
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Using the expression of the rates given in Eq. we have d^w^m, m ± 1) = 

(0(??2 ± 1) — 0(m))P/2 ~ ±f3ed x O(x)/2. Inserting this in the above equation and using 
Eq. IA.81 one obtains the following expression of the local current v t < in the continuous 
limit 

W ( X ) = EiM dx0 ( x ). (ah) 

Pt'ix) 

Together with the first term in the response given in Eq. HSJ one recovers from this the 
response function given in Ref. [33], namely 



R(t, if) = P^- (A(x t )0(x t ,)) - ( A {x tl )^p\d Xt ,0{x tl 
at' \ Pt>{xt') 



(A.12) 



Appendix B. Generalization to multidimensional linear Langevin dynamics 

We now extend the results obtained in section 13.11 for Langevin dynamics with one 
degree of freedom to many dimensions. To that end, we start with the following 
multidimensional linear Langevin equation: 

x t = N t x t + h t + rj t , (B.l) 

where r\ t is a white noise such that (rj t ) = and (r)tr)t>) = 2T t T t S(t — t'). We denote T t a 
time-dependent temperature, T t a positive and symmetric matrix, and N t an arbitrary 
matrix. 

For any t > s, we introduce the following matrix 

T(t, s) = exp ( f duN u ) , (B.2) 



where the exponential is to be understood as time-ordered. This matrix satisfies the 
useful identity that for any t > s > u, T(t, s)T(s,u) = T(t,u). We also introduce the 
matrix 

D(t,s)=2 I duT{t,u)T u T u T\t,u), (B.3) 

J s 

where the dagger denotes the transpose of a matrix. 
The solution of Eq. IB. II is 

x t = T{t,0)x + [ dsT{t,s){h s + r] s ). (B.4) 







We assume that the initial condition xq is distributed according to a Gaussian, 
with a characteristic mean m = (xo), and an initial covariance matrix Vq = 
(xq — (xo)T (xq — (xo)y ). From Eq. IB. 41 we deduce the mean 



m t = (xt) = T(t, 0)m + / dsT(t,s)h s , (B.5) 

Jo 

and the symmetric matrix of the covariance 

V t = ((x t - mt) ix t - mt)) = T(t, 0)V o T*{t, 0) + D(t, 0). (B.6) 
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As in the ID case, we note that in this case too, V t is independent of the perturbation 
h t . Since x t is a Gaussian random variable, its distribution is 

P t (x, [ht]) = (det (2ttT4))-5 exp (-~(x - mjVf^x - m t )J , (B.7) 



and therefore 



n t (x, h t ) = (det (2tt^)H exp ( - \{x - m t )V t ~\x - m t ) ) , (B.8) 



with 

m t = T(t : 0)m + ( f dsT(t,s))h t . (B.9) 



Since 

d ht {lmr t {x, h t )) 



(J^ dsT\t, s)) Vf\x - m t ), (B.10) 



The functional y of interest here, takes the form 



so that 



yr = -J o dth < [I dsT] (^ s )) V f\ x t ~ ™t)- (B.ll) 
From this, we obtain the average 

(y T ) = - dth tQ o dsT^t,s)^V t - l (m t -m t ), (B.12) 
= jf dth t ^jf dsT\t,s)^j Vf 1 ^jf duT(t,u) (h t - h u )^j , (B.13) 

y T _ (y T ) = _ jf ^ ^ jf dsT t( t> y r i( x< _ mt ). (b.14) 

Through explicit calculation, one can verify that 

(OV - 0V» CVr - OV»> (B.15) 

= 2^ <M 4 ^jf dsT\t,s)^j Vf 1 J q dt'T(t,t') ^jf dsT(t',s)J 
Using an integration by parts, one can check that 

jf dt'T(t, t') ^ dsT(f, s) j ^ = jf dt'T(t, f) {h t - h t >) . (B.16) 
Taken together, these two equations imply the relation 

({y T - (y T )) (y T - (y t ))) = 2 (y T ) . (B.17) 

Since y? is a Gaussian variable, the detailed fluctuation relation given in Eq. 1661 follows 
from this for the general linear Langevin dynamics. 
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